Spin connection and boundary states in a topological insulator 
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We study the surface resistivity of a three-dimensional topological insulator when the boundaries 
exhibit a non trivial curvature. We obtain an analytical solution for a spherical topological insulator, 
and we show that a non trivial quantum spin connection emerges from the three dimensional band 
structure. We analyze the effect of the spin connection on the scattering by a bump on a flat 
surface. Quantum effects induced by the geometry lead to resonances when the electron wavelength 
is comparable to the size of the bump. 
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I. INTRODUCTION. 

Strong Topological Insulators (TI 's) are a new class 
of materials with a bulk gap but surface states denned 
on surfaces of all orientations 1 - - — , making the boundaries 
gaplcss. The number of surface states at a flat surface 
with a given orientation is odd, and each of them shows 
a conical singularity, described by the two dimensional 
Dirac equation 4 -. Localized states also exist at other lat- 
tice defects, such as screw dislocations*—. 

The transport features of electrons at the surfaces of 
TI 's is being intensively studied. The wavefunctions 
have an internal spinorial structure made up of two 
slowly varying components related by time reversal in- 
variancc. Backscattering due to smooth perturbations 
which preserve the time reversal symmetry is forbidden, 
making the transport properties of these compounds sim- 
ilar to those of graphene in the absence of intervalley 
scattering*. Surfaces with a finite curvature allow for 
scattering processes due to the existence of a non triv- 
ial metric, which has been studied in the classical limit 9 , 
when wavepackets are well approximated by point parti- 
cles following classical trajectories. 

An analysis of the electronic properties of curved sur- 
faces of TI 's requires information about the way a non 
trivial metric changes the effective Dirac equation. It 
is well known that the Dirac equation could be written 
on a curved space time introducing the spin connection 
and the rotation of Dirac matrices* -. The existence of 
the spin connection has been postulated in topological 
insulators*^. The emergence of the spin connection from 
the combination of the three dimensional electronic struc- 
ture of a TI and the two dimensional metric of a boundary 
with intrinsic curvature has not been studied so far. 

In the next section, we analyze the surface states for 
the simplest curved boundary with a non trivial metric, 
the sphere. The conservation of the angular momentum 



in this geometry allows us to calculate the entire spec- 
trum of surface states, and to show that the spin connec- 
tion term is induced in the effective surface hamiltonian. 
We use this information in Section III to analyze the ef- 
fect of the curvature in the scattering by a bump in a flat 
surface, a process considered in the classical limit in 9 -. 
Related processes can be defined in graphene with topo- 
logical defects* 2 -. Technical details of the calculations are 
discussed in the appendices, including an analytical study 
of the boundary states in a cylinder, calculated numeri- 
cally inia. The boundary of a cylinder can be considered 
a surface without intrinsic curvature and spin connection. 



II. MODEL OF A SPHERICAL TOPOLOGICAL 
INSULATOR. 



The surface states of model single particle hamilto- 
nians for a TI have been studied particularly for a 
flat boundary*^ and for an infinite cylinder boundary 
surfaces*. A minimal model reproducing the band struc- 
ture of a TI requires four orbitals, related in pairs by the 
time reversal symmetry**. A further simplification in- 
cludes just the linear in-momentum contributions to the 
hamiltonian of Rcf. Il5f 



H = 7°A + hvpfh 



(III) 



where vf is the Fermi velocity and the matrices j a 
are given in terms of Pauli matrices by 7 = tq ® t z , 
7 1 = cr x ® t x , 7 2 = -<7 y ®t x 7 3 = a z <E> t x . Here a a and 
Tb denote matrices in the spin and even-odd orbital par- 
ity spaces, respectively. This hamiltonian satisfies time 
reversal symmetry T = K,ia y <8>l2x2 (here I2X2 the 2x2 
identity and K, the complex conjugation). Bulk eigen- 
functions in cartesian coordinates are: 
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(k± — k x ± ifcy), where the band energies are e±(k) = 



1,2 

K y 



fcf ) and iV± is the norm of the 



'A 2 
states. 

Surface states appear in this model if the gap param- 
eter, A, changes its sign at the boundary, so that, e.g., 
A > in the inside, and A < in the vacuum. The 
model allows also for the analytical computation of the 
surface bands of a cylinder—, as shown in Appendix [A] 
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Later we will include also quadratic corrections to the 
hamiltonian in Eq lII.ll and we will show that boundary 
conditions need to be chosen in a different way in that 
case. 



In order to obtain the solution of the hamiltonian III, ll 
onto a sphere we rephrase its eigenvalue equations into 
spherical coordinates. The eigenvector of energy E, ^ = 
{^fA, *c, *d), satisfies the equations: 
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(II.3) 



(here H = vf = 1), with the boundary conditions 

A ^) = {a1 III ( IL4 ) 

We choose A^„ = —A out = A for simplicity, so that 
the exponential decay of the boundary states into the 
bulk near a flat surface is defined by the length scale 
A = hvp/A. The angular momentum is conserved and is 
quantized in half integer units (see, for instance Ref. [ill ). 
Its eigenfunctions allow us to reduce the set of Eq.s pi.3p 
to two coupled differential equations for the radial co- 
ordinates, as discussed in Appendix B. It can be shown 
that the energy spectrum converges exponentially to the 
one of the two-dimensional Dirac equation onto a sphere: 



Ej = ± 



hv F (J+l/2) 
R 



1 + O e- R / A 
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(II.5) 



where J m ax ~ R/A. The multiplicity of each level is 
2 J + 1. The exponential convergence of the energy levels 
to the asymptotic value in eq. III.5l is shown in Fig. [T] This 
type of convergence implies that the effective hamiltonian 
describing the surface modes does not admit an expan- 
sion on higher order derivatives, of the type A(Adi) n . 



The study of the hamiltonian in cq. III. II can be extended 
in a straightforward way to the case when A out ^ Aj n , 
although it becomes cumbersome to obtain analytical 
expansions in the limit R — > oo. The numerical solu- 
tion, obtained by generalizing the analysis of Appendix 
B, shows an agreement with the spectrum in Eq. III. 51 of 
the same accuracy as those reported in Fig. [1] 

As quadratic terms do not break the spherical symme- 
try, they can be safely added to Eq. III. II by the simple 
substitution A —> A + a(k\ + ky + k^), where a is a 
constant. Hence the angular part of the wavefunctions 
remains unchanged while its radial part satisfies second 
order coupled equations, in place of those in Eq. III. 31 
For each value of the energy, E, we find evanescent 
waves with two different decay lengths, Ai(£7) and A2(E), 
which are given by the roots of a fourth order polyno- 
mial. The boundary conditions need to be replaced. The 
simplest boundary condition compatible with the new 
second order equations is ^a(-R) = ^b(R) = ^c(R) = 
^d(R) = Ol^. By solving numerically these boundary 
conditions, we find again an agreement with eq. |IL5] sim- 
ilar to that shown in Fig. [TJ Results are shown in Fig. [3] 

We conclude that the boundary states on a spheri- 
cal TI satisfy the Dirac equation on the surface of the 
sphere. The spin connection, related to the intrinsic cur- 



3 




Log[\Ej\x(R/(J+l/2))-l] 
-3 
-4 
-5 

2 4 6 iS 10 
Log[\Ej\x(R/(J+l/2))-l\ 

-5.0 

-5.5 

-6.0 

-6.5 



to a local incrtial frame: 



1 2 3 4 5 6 7 



J 



FIG. 1: (Color online). Dependence of the surface energy 
levels on angular momentum, J, and radius, R. The deviation 
from the result for the two dimensional Dirac equation on a 
sphere, see eg. III. 51 is plotted. Top: Dependence on R. From 
top to bottom, J = 1,2,3. Bottom: Dependence on J. From 
top to bottom, R = 8, 10, 12. In all cases, vf = 1 and A = 1 
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FIG. 2: (Color online). Energy levels of a spherical topologi- 
cal insulator of radius R = 30 with a quadratic dispersion re- 
lation, obtained by the replacement A — > A+a(k^.+ky+ki) in 
eg. III. II Other parameters are A = 1 and Vf = 1. The bound- 
ary conditions are ^a(R) = ^b{R) = ^c(R) = *£>(#) = 0. 



vature of the metric, clearly emerges at the boundaries 
of a TI. More generally, the boundary states satisfy the 
Dirac equation on a curved space-time^ 



7 "(0„ + r /1 )* = o 



Here 7 M 

7 Q = 2 = - 



l j a ( with 7 
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(II.6) 

a=1 = cjV and 



-o x ) are the rotated Dirac matrices, satisfying 



the generalization of the flat algebra 



(II.7) 



The Minkowski metric rj a b has been replaced by the 
curved one g^ v — rj a b e a e b v . The tetrads e a are the 
elements of the Jacobian matrix of the transformation 
from the coordinates x M , defined on the whole manifold, 



dx a 



(II. 



|r% b E ab where E a6 



is the spin connection 
«/2[7 ,7(,] are the generators of the spinorial representa- 
tion of the Lorentz group, and the connection coefficients 
T a b are given by T a b = 



III. SCATTERING OFF A GAUSSIAN BUMP 

A. Unrelaxed lattice 

We now derive the resistivity for electrons propagating 
at the fiat boundary surface of a TI, when they are scat- 
tered off a gaussian bump of height z(\r\) = h e~ r . 
The Boltzmann relaxation-time approximation can be 
used (f(0) = kp I (irhvp) is the density of states at the 
Fermi level for both spins): 



p{k F ) 



1 



e 2 v 2 F is(0) r(kp) 



(111.9) 



where the usual definition of the total relaxation rate is 



2- 



r(k F ) 



i^(o) [ de(i-k-k'^\(k\t eif \k')\ 2 {iii.iQ) 



Here (k\t e ^\k 1 ') is the matrix element of the t— matrix, 
which depends on the energy and on the scattering angle 
between the incoming and outgoing wave. 

Since the metric induced on the manifold by the bump 
is axially symmetric, it is convenient to rewrite the two 
dimensional Dirac equation in flat space time in cylindri- 
cal coordinates: 



ihvp [ a r d r + a 6 -de ) * = E$ . 

r 



Here the matrices a r ' are a r 
a 6 = — sin# cj x + cos9 a v . 
Given the metric^ 



-10 
g^=\ 1 + f(r) 



(III.ll) 
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(111.12) 











where f(r) = (dz(r)/dr) 2 , we rewrite Eq. (jII.6[) and, 
as shown in Appendix C, the spin connection 
can be embodied in the wavefunction as a real pref- 
actor \& = $ exp J r + °° dr' Ag(r'), where A$(r) = 

The spinor $ satisfies the equation 



(yr+jld-i). 



d 



,1 d 
er -■ 



$ = sk $ 



(111.13) 



y/l + f(r) dr r 89 
where s = +(— ) for particles (holes) and E — hvpk. The 
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real prefactor can be interpreted as the origin of charge 
puddles accumulating at the bump. Eq. 13[) de- 

scribes an unrelaxed lattice. Relaxation of the structure, 
besides adding an effective gauge potential, may further 
change the spin connection. As elastic deformations do 
not add any curvature, the change only implies a trivial 
holonomy on the wave function. This is a way of restat- 
ing the Saint Venant conditions for the two-dimensional 
case. Changes in the Dirac Equations are well localized 
in space close to the bump, hence a scattering picture 
can be fruitfully adopted here. We will focus on the par- 



ticle sector of the theory, and assume that the incoming 
A:— vector is in the direction of the polar axis (6 = 0). 
The eigenfunctions can be expressed as superposition of 
angular momentum m eigenstates: 



Mr,0) = 




(III. 14) 



The Born approximation is worked out, to lowest order, 
in Appendix C. Using an asymptotic expansion of the 
wavefunctions given in Eq. (|D.9[) . we have: 



(k F ) 



k F X ^ 



[sin 2 8 m - cos (5 m+ i - 5 m -i) sin<5 m+ i sin(5 m+1 ] (III. 15) 



where the phase shifts S m for the m th component of the 
wavefunction are reported in appendix C and rib is the 
concentration of bumps. 

At low incoming electron energy, it turns out that the 
terms with m — 0, ±1, ±2 are 0[(k£) 4 } and when choos- 
ing 4TT 2 h 2 /£ 2 ~ 1, they sum up to S « 0.733. The terms 
with m = ±3 are 0[(k£) 8 ], while the terms m = ±4 are 
0[{k£) 12 } (see Appendix C). 

Eventually, the resistivity for independent point like 
defects, when the carrier density is low, (i.e. low incom- 
ing energy) is: 



All other terms scale for m >> ki as 



p{k F ) 



^n b n£ 2 is(j) (k F £) 2 + 0[(k F £)^ 



(111.16) 

S(h/£) is a numerical prefactor which depends on the 
strength of the perturbation parametrized by h/£. The 
plot of the resistivity vs energy in dimensionless units, 
k F £ 1 for various values of h/£ is shown in Fig. [31 The 
leading term is proportional to the density of carriers n. 

It has been proven recently that the classical limit for 
large incoming energy (i.e. relatively high densities n) 
corresponds to an energy independent v F r(k F )Q. This 
implies that, according to Eq. pil.9|) . p oc h/(e 2 k F ) ~ 
l/y/n, in this limit. We derive the same conclusion with 
a careful analysis of the sum in Eq. (|III.15[) at k£ >> 
1. Classically, angular momentum conservation in the 
scattering implies that m ~ kb. Here b is the impact 
parameter measured from the center of the bump in the 
direction orthogonal to k. By asymptotically expanding 
the Bessel functions appearing in the phase shifts, it turns 
out that there is a collection of terms contributing to the 
sum, which are roughly independent of m, as long as 
k£ » m. For these m values, tan5 m is of the form: 



tan 5„ 



"7 



M 
2 



2 1 ' {k£f 



(111.17) 
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(111.18) 



and therefore they rapidly converge to zero. 

We conclude that in the semiclassical limit W > 1 a 
factor k F comes from the relevant terms in the sum of 
Eq. piI.15[) . which are all of the same order. This fac- 
tor cancels with the k F appearing in the denominator, so 
that the result for v F r(k F ) is independent of k F . This 
is in fact found numerically. In Fig. |3j we see that the 
conductivity ~ p^ 1 grows linearly at large k F £ and has 
a minimum in the neighborhood of k F £ ~ 1. In Ap- 
pendix D we report a simple argument based on a saddle 
point approximation of the ^2 m which qualitatively re- 
covers the classical limiting result for large k F £, derived 
in Ref. |9|. In Fig. [3] it is shown a significant increase of 
the cross section of the bump for k F £ ~ 1 and h/£> 0.2. 
This increase is due to quantum resonances induced by 
the non trivial spin connection. 



IV. DISCUSSION. 

We have shown in Section II how the two dimensional 
Dirac equation in curved space emerges at the simplest 
boundary with non trivial metric, the surface of a sphere. 
The metric enters through the spin connection, which re- 
flects the properties of the internal spin under parallel 
transport along the surface^. The spin connection re- 
flects a quantum feature of the electrons, and cannot be 
inferred from solely classical arguments. The spin con- 
nection leads to a finite Berry phase when the electron 
is transported around a closed geodesic. A manifesta- 
tion of this effect is the quantization of the total angular 
momentum in half integer units. 

The model that we have studied leads to simple ana- 
lytical expressions of the energies and wavefunctions of 
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We have analyzed the scattering of Dirac fermions by 
surface corrugations which induce a non trivial curvature 
in the quantum limit, kpi < 1. We find that the resistiv- 
ity due to a finite concentration of bumps, rib, vanishes 
as kp for small fejr, due to a combination of a density 
of states factor, which goes as kp, and a scattering time 
which increases as fc^ 3 . By comparison, the scattering 
time in the classical regime^ — > oo) is independent 
of kp, and p ~ kp . The wave nature of the quasipar- 
ticles allow them to diffract around the bump, making 
it effectively transparent for long wavelengths, kp£ <C 1. 
The non trivial curvature induces quantum reseonances 
for kpl ~ 1 and an aspect ratio h/i > 0.2. 



FIG. 3: (Color on-line) Main Panel: inverse scattering time 
as a function of kpl for different values of the aspect ratio 
of the bump h/l. Inset: the resistivity due to scattering off 
a bump on the surface of a topological insulator, in units 
p-e 2 A/(2ivM 2 ) vs. dimensionless energy k£, at different ratios 
h/i. 



the boundary states. They can be used as a zeroth ap- 
proximation to situations close to spherical symmetry, or 
where an isotropic electronic structure can be obtained 
by rescaling a length. We find that the corrections to the 
two dimensional Dirac equation depend exponentially on 
R/A = RA/HvpH. 

The model describes external surfaces of mesoscopic 
crystals and internal voids in bulk systems. In the case 
of a small void, we find that two doublets at a finite 
distance of the Dirac energy appear for radii R > 2/ A. 
These voids will act as molecules embedded into the 
bulk material. The interaction energy of electrons lo- 
calized inside the voids scales as the level separation, 
Eim ~ e 2 /(ei?), where e is the dielectric constant of 
the topological insulator. At temperatures below this 
scale, voids with an odd number of electrons will give 
rise to magnetic moments. The RKKY interaction be- 
tween moments at different vacancies should decay ex- 



ponentially, J rkky {y - ?' 



-|r r |/A_ Hence, small 



vacancies might give rise to a paramagnetic susceptibil- 
ity in topological insulators. If they are within a dis- 
tance d ~ A from the surface, these local moments will 
hybridize with the surface states, leading to the Kondo 
effect^. 
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Appendix A: Boundary states at cylinder surface 

In the following appendixes we will use h = vp = 1, 
except in main results. Let us start from the k ■ p model 
Hamiltonian of Eq. (III.ip : 



H[r] 








id z 

-A i(d : 
i(d x - idy) A 
-id 



+ id v 



i(d x + id v )\ 


-id, 
-A 
(A.l) 

To find surface states in this approximation is enough 
to match the solutions of the Schrodinger equation at 
the surface of the cylinder. The gap A should change 
its sign between in and out of the surface. We rewrite 
the eigenvalue problem in cylindrical coordinates for the 
4-component spinor ("F^, >Fb, ^c, ^d)' 
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(25 - A) ^( ? ) = i|* B (?) + e " - 

(25 + A) * B (r) = i|^(r) + e» - ~) *c(F) 

(25 - A) * c (r) = e"« (i£ + ^) * B (f) - i|^(f) 

(25 + A) * D (F) = e" ie ^ + ~\ 9 A {f) - i£*c(?) (A.2) 

(in the following k is in the z direction, which is the axis of the infinite cylinder. Inside the cylinder, the wavefunctions 
which arc mostly localized close to the surface involve the modified Bessel functions I n {nr) with integer n. They diverge 
exponentially at infinity but are finite for r — > 0. The two eigenfunctions at fixed energy 25 are (k is unknown for the 
moment): 

/(25 + A) /„(/«■) \ / iKl n (Kr)e ie 

l^'^-JV e e ,|25,2<)-- _ fc/n+l(Kr) |e { >e . (A.3) 

\ lK 2 n+1 («r) e- 4 7 V(25 - A) J n+1 («r) 

The energies of these states are 25 = ± \J A 2 + k 2 — k 2 . Outside the cylinder, the functions K n (nr) replace the I n (nr), 
as the former decay exponentially for nr — > oo and A — > — A. The eigenfunctions are: 



(25- A) K n (nr) \ / -i« 2sT„(Kr) e ie 

|p 1 \ _ 1 k K n (hir) „-in6 -ikz P n i 1 I -i(„+l)e -ikz 1 a j-v 

l^' 1 ^-/?! e e ,|25,2>)-- _fc^ n+1 ( Kr ) l e ' ' e ■ ( A - 4 ) 
,- iK 2C n+1 ( K r)e- 4 V \(25 + A)2f„ +1 ( K r) 



The eigenvalues are again those of Eq. ([A.3[) matching conditions at R, the radius of the cylinder, lead 

The two wavefunctions inside the cylinder should be to 
matched to the two outside for each value of n. The 



Dct 



inl n (kR) (E + A)I n (nR) 

kI n (nR) 
-kI n+1 {KR) 



-iKK n {nR) (25 - A)K n (nR) 

kK n (nR) 
-kK n+1 (nR) 



(25- A)I n+1 ( K R) iKl n+1 (nR) (E + A)K n+1 {nR) -iKK n+1 ( K R) 
The vanishing of the determinant implies: 

[l 2 n K 2 n+1 + Kl I 2 n+1 ] k 2 + (2k 2 - 4A 2 ) I n I n+1 K n K n+1 = . 

I 



(A.5) 



(A.6) 



The presence of products of the Bessel functions 2„ and 
K n assures that, in the limit of kR » 1, there is just 
an inverse powerlaw dependence of the secular problem 
on kR. and Eq. (|A.6p becomes: 



Hence, the energy of the states is, according to Eq. (|A.8[) 



hvf 



(A.8) 



= A 



1 - 



(n + l/2) s 
A 2 2? 2 



(A.7) 



This result is in complete agreement with Eq.5 of Rcf. 
[l3l . Let us now consider the opposite limit kR « 1. 
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Expansion gives, up to second order in 1/nR: 



The energy reads in this limit: 



k 2 \ 2n(n+ 1) 



16 n 2 (n + iy 



A 2 . (A.9) 



E n {k) » ±A 1 



1/ 2n(ra+ 1) 



16 n 2 (n + l) 2 



/ hv F k 2 



r 



Appendix B: Angular momentum eigenstates 

Generalized angular momentum operators J can be 
defined as usual as the sum of spin and orbital angu- 
lar momentum. It can be shown that the Hamiltonian 
in eq. III. 31 commutes with J 2 , J z therefore its eigenstates 
can be labeled by \j,m) with J 2 |j, m) = j(j + 1)\J,J Z ) 



J 



and J z \j,m) = m\J, J z ). In order to obtain single valued 
eigenfunctions, the values of J and J z must be half in- 
tegers. As usual by using J + \j]j) = and J~\j,m) oc 
\j, m — 1) we can explicitly construct the wavefunctions 
of the different states \ j, m) . The Hamiltonian eigenfunc- 
tion can be thus expanded onto the the lowest angular 
momenta states: 
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\ sin 2 (0)e 2 ^ 
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y 3sin(0)cos(0)e^ J 



(B.l) 



sin(0)e 



where the states 1 3/2, — 1/2) , 1 3/2 , — 3/2) are not explic- 
itly exhibited, here. 



form: 



Appendix C: Spherical boundary states and energy 
spectrum 



It can be snown that boundary states in the spherical 
case for j = m = n— 1/2 (n > 0) have the following 
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/i T (r) 
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/sin"- 1 ee^ n -^\ 







(C.l) 



Here f~(r) and f + {r) are radial functions localized at 
the boundary for r < R and r > R, respectively, and 
they satisfy the equations: 



(E T A)/+ = -idrf? --(n + 1)/+ 



(E±A)ff = -id r f? + -( n -l)f? 



(C.2) 



The system can be decoupled in a pair of Bessel equations 

n(n+l)~ 



dr 2 1 r dr 1 
dr 2 2 r dr 2 



(A - E 2 
(A 2 - E 2 ) 



n(n — 1) 



fx = , 



ft = o , 



(C.3) 



whose solutions are: 



/f(r) 



• A — E 



C i n (nr) 



(C.4) 



z/ r < i? 
if r > R 

where i n , fc„ are the modified spherical Bessel functions: 



J2 ( r ) = C in-i(nr) 

' f+(r)=-iC+ ^k n {Kr) 
f+(r)=C+k n _ 1 ( K r) 



in(x) = y^4+i (x) , fe„(x) = <J^-K n+ i(x). (C.5) 
The matching conditions can be written written using 

E q . jcaji 



-iC~ (A - £7) £„(«#) = -£C+ (A + E) k n {nR) 
C~ i n -x(KR) = C+ fc„_i(/c.R), 



which give rise to an implicit equation for the eigenener- 
gies of the system: 



A - E _ k n {nR) i n -x(KR) 
A + E ~ i n {nR) k n -i(nR) 



(C.7) 



This equation, in the limit AR — > 00, gives the admissible 
values of the energy: 



En = ^ n ~lt ' n = 1, ... ,n max (C.8) 



which is reported in Sec. II. 



Appendix D: Derivation of the elastic (-matrix for 
scattering off a gaussian bump 



We derive the elastic t— matrix when a localized de- 
formation is present on the surface of a TI. As ir>i£ we 
consider a two-dimensional spatial sheet modeled on a 
two-dimensional axial symmetric manifold with a single 
gaussian bump. The axial symmetric gaussian surface 
may be represented in Minkowski space-time by the func- 
tion 4> = (t> x, y, h(r)) with r 2 = x 2 + y 2 . From Eqs. 
pil,12[) . we may read the tetrads 



sin 9 



2 sin# 2 cos# 



(D.l) 



(C.6) The Dirac equation on a radially symmetric manifld is 



d r + a e \ -d e 



2r 



# = E^ 



(D.2) 



r 



The gauge potential in (|D.2|) is the spin connection 

1 



We pose ^ = <& exp f dr' Ag(r') with Ag is the spin 
connection above. The m component of the spinor $ has 



r *=2 1 



<J Z 6 



(D.3) 
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the form (m is the angular momentum integer): 



<f> m (r,6\k,s) 



u sm {r) 
is v sm (r) e v 



(D.4) 



ing (|D.4[) in the Dirac eq. pil.131) and dropping the labels 
sm, we find that the functions u(r),v(r) have to satisfy 
the following equations: 



where Ok is the angle that the direction of the k vector of 
the incoming wave forms with the polar axis. Substitut- 



J 



1 d 2 u{r) 1 du{r) 



y/1 + f dr 2 r dr 
1 d 2 v(r) 1 dv(r) 



1 



dr^TTJ 



du(r) 
dr 



■y/T+fu(r) + k 2 u(r) 



y/l + f dr 2 r dr 



1 



dv(r) (m + 1)^ 



dr 



y/T+f v(r) + k 2 v(r) 



0. 

= 0. 



(D.5) 



Due to the symmetry of the problem is suitable to expand 
the Green's function for the flat space-time problem in 
polar coordinates 



2- 



e™ (e - 9 'Wr,r'). (D.6) 



be the gaussian bump z(|r|) defined in sec. III. A. This 
implies that /(r) = (4h 2 r 2 /£ 4 )e~ 2 ( r /^~ . We assume that 
the ratio h/i is small, so that we can expand Eq.s (|D.5|) 
by retaining just the lowest power oih/l. By comparison 
with the system for the flat space (i.e. f(r) = 0), we 
define the perturbative potential : 



The Green function displaying the correct jump of the 
derivative at r = r' is: 



g m (x,x') = 2ir 2 J m (;£<) Y m (x>) 



(D.7) 



Here r< (r> ) is the smaller (larger) of the two arguments 
r,r'. We now specialize the shape of the bump h(r) to 

I 



dr 2 



4 8r 
r~¥ J dr 



(D. 



In the Born approximation, the Dyson equation for 
e.g. Ukm reads: 



2 roo 



Ukm(r) = J m (kr) + —J dr' r' g m (r, r')V m (r') J m (kr') 



We have used the fact that J m solves the Bessel differ- By defining : 
cntial equation to simplify the action of V m on J m itself. 

I 



(D.9) 



2I.2L2 r°° 



tan 5 ri 



A>K 2 k z h 



r 



(D.IO) 



the scattering state for r/t — > oo takes the form Mfe m (r) ~ u v, m •<-> — m— 1, it is easy to see, that the sums which 
Jm(kr) +tan<5 m Y m (kr). By exploiting the symmetry of include S m , for all m, arc equal. Therefore, our result is 
the Dirac masslcss equation with respect to replacements valid for both components of the spinor solution given by 
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Eq. piI.14[) . implies that, far from the bump ( rjl — > oo ), the outgo- 

The integrals of Eq. p3.10p can be evaluated analyti- ing wave takes the form: 
cally. The asymptotic expansion of the Bessel functions 



1 tan<5 m . . ran ir 

Ukm{r) ~r«-S-oo ; = = COSXm H , „ = Sin Xm = COS (Xm + Omj Xrn = kr - 

V/l + tan 2 6 m Vl+tan 2 (5 m 2 4 



I 

We now evaluate the t— matrix element for a scattering scattered elastically by the gaussian bump and a plane 
event, in which an incoming wave with wavevector k is wave of wavevector p emerges. The t-matrix element is: 



(p\t(k)\k) 



TkR 2 



i7r/4 



R 



dr 



■ f r R dr' r' 2 



3 -im(e p -e fc ) 



The space integral arises from the exponential prefactor 
of \I/ defined before Eq. pil.13p . To evaluate the relax- 
ation time formula of Eq. piI.15|) we first perform the 
integral over the angle 9 P of the square modulus of the 
angle dependent exponentials in the sums. The integral 
is non vanishing only for m — ml = ±2, 0. By rearranging 
the sums then eq. (|III.15[) is obtained. 



Appendix E: Semiclassical approximation 

In this Section we present a tasteful derivation of 
the classical high energy limit for the relaxation time. 
The latter can be obtained by assuming classical diffu- 
sion along the geodesic trajectories across the bump and 
yields^: 



Here b is the impact parameter of the incoming particle, 
while 9 is the scattering angle and A is the area of the 
sample. The starting point is the usual expression for the 
relaxation time of Eq. piI.15[) : 



kA 



d9 (1 



(E.2) 



given in terms of the scattering amplitudes f m (6) = 
r e i25 TO _ -n e ime _ hjgh energy, many m terms con- 
tribute to the sum, so that we take its continuum limit, 
which amounts to integrate over continuous values of the 
classical angular momentum m = kb. As forward scatter- 
ing is excluded form Eq. pi!.2[) . it is enough to apply the 



saddle point approximation to the resulting integral 2 ^: 

^ e i(2* m+m 0) ^ e i(26 mo+ ra 9) x I dme ^\j m - m ^\ 

(E.3) 

Here mo is the stationary point, which solves the saddle 
point equation: 



d8 m 
dm 



0. 



(E.4) 



Derivation of this equation, once more, provides a rela- 
tion between the second derivative of the phase shift 5 m 
and the the angle 9 



d 2 S r , 



1 d9 



dm 2 2 dm 



= 0. 



(E.5) 



An analytical continuation in the complex m plane allows 
us to make the integral converge. Gaussian integration 
in Eq. pil.3p implies that 



dm 



dO 



(E.6) 



thus yielding the expected result: 



1 VF_ 1_ 

t k A 

ta r m y 



J d9(l - cos( 



dm 



dO 



VF 

A 



db 9 2 {b) . 



(E.7) 



In the last equality the conservation of the angular mo- 
mentum m = kb has been exploited, together with the 
remark that the scattering angle only depends on b in 
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the classical diffusion. This reproduces the desired high so that we recover again 
energy behavior. 

The analysis of Eq. (|III.15| provides a similar conclu- 
sion. The quantity <5 m +i — <$ m -i, appearing as the argu- 
ment of the cosine, is w 2 d5 rn /dm. At large incoming 

energies, 6 m « (2m + l)n/2, which is consistent with the \ Vp r (p- Vp r 

asymptotic form of the wavefunction given in Eq. (|D.lll) . — ~ / ^ m ~^ ~ J ^ ■ (E.9) 



To lowest order we get, according to Eq. (|E.4 

sin 2 5 m - cos (S m+ i - 5 m -i) sin<5 m+ i sinJ m _! « — , 

(E. 
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